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Abstract 

It's well-known that there is a very powerful error bound for Gaussians 
put forward by Madych and Nelson in 1992. It's of the form|/(x) — s(x)\ < 
(Cd)d \\f\\ h where C, c are constants, h is the Gaussian function, s is the in- 
terpolating function, and d is called fill distance which, roughly speaking, 
measures the spacing of the points at which interpolation occurs. This 
error bound gets small very fast as d — > 0. The constants C and c are 
very sensitive. A slight change of them will result in a huge change of the 
error bound. The number c can be calculated as shown in [9] . However, 
C cannot be calculated, or even approximated. This is a famous question 
in the theory of radial basis functions. The purpose of this paper is to 
answer this question. 
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INTRODUCTION 

Let h be a continuous function on R n which is conditionally positive 
definite of order m. Given data {xj i fj),j = 1,...,N, where X = {xi ...,xjf} 
is a subset of points in R n and the /js are real or complex numbers, the 
so-called h spline interpolant of these data is the function s defined by 

N 

S(x)=p(x)+ £ Cjh(x-Xj), (1) 

where p(x) is a Polynomial in V m -\ and c ' s are chosen so that 



N 

E Cjq(xj) = (2) 

3=1 
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for all polynomials q in V m -\ and 



N 



P( x i) + E Cjh(xi - xj) = fi, i 
3=1 



= 1,...N. 



(3) 



Here V m ~\ denotes the class of those polynomials of R n of degree 

< m — 1. 

It is well known that the system of equations (2) and (3) has a unique 
solution when X is a determining set for V m -\ and h is strictly condi- 
tionally positive definite. For more details see [7] . Thus, in this case, the 
interpolant s(x) is well defined. 

We remind the reader that X is said to be a determining set for V m -i 
if p is in T m —\ and p vanishes on X implies that p is identically zero. 

1.1 A Bound for Multivariate Polynomials 

A key ingredient in the development of our estimates is the following 
lemma which gives a bound on the size of a polynimial on a cube in R n in 
terms of its values on a discrete subset which is scattered in a sufficiently 
uniform manner. For its proof, please see [9]. 

LEMMA1. For n = 1, 2, ...,define 7 n by the formulas j 1 = 2 and, if 
n > 1, 7 n = 2n(l +7 n _ 1 ). Let Q be a cube in R n that is subdivided into q n 
identical subcubes. Let Y be a set of q n points obtained by selecting a 
point from each of those subcubes. If q > j n (k + 1), then for all p in Vk 



1.2 A Variational Framework for Interpolation 

The precise statement of our estimate concerning h splines requires a 
certain amount of technical notation and terminology which is identical 
to that used in [8] . For the convenience of the reader we recall several 
basic notions. 



sup|p(x)|< e 2n7n ( fc+1 )sup \p(y)\ ■ 
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The space of complex valued functions on R n that are compactly sup- 
ported and infinitely differentiable is denoted by D.The Fourier transform 
of a function (f> in V is 

= J e- i<x ^>4>(x)dx. 

In what follows h will always denote a continuous conditionally posi- 
tive definite function of order m. The Fourier transform of such functions 
uniquely determines a positive Borel measure /j, on R n ~ {0} and constants 
a 7 , |7| = 2m as follows: For all if) G V 

fh(x)^(x)dx= f\h)-x(0 £ Ww^ldMO 

J [ \-y\<2m J 

+ £ DT$(0)if, (4) 

| 7 |<2m 

where for every choice of complex numbers c Q , |a| = m, 

E £ a a+ pc a cp > 0. 

|a|=m|/3|=m 

Here x is a function in V such that 1 — has a zreo of order 2m + 1 
at £ = 0; both of the integrals 

/o^KiKI*"^), feidMK) 

are finite. The choice of x affects the value of the coefficients a 7 for 

I7I < 2m. 

Our variational framework for interpolation is supplied by a space we 
denote by C h , m - If 

T^m = {<t> € V : f x a <j){x)dx = for all \a\ < m} 
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,then Ch t m is the class of those continuous functions / which satisfy 

\f f(x)<i>(x)dx\ <c(f){fh(x-y)<l>(xWyjdxdy} h (5) 

for some constant c(f) and all (f> in V m . If / € Ch, m let \\f\\ h denote the 
smallest constant c(/) for which (4) is true. Recall that is a semi- 

norm and Ch jTn is a semi-Hilbert space; in the case m = it is a norm and 
a Hilbert space respectively. 

Given a function / in Ch, m and a subset X of R n , there is an element 
s of minimal C\^ m norm which is equal to / on X. 

The function space Ch, m is not very easy to understand. However Luh 
made a lucid characterization for this sapce in [3] and [4]. For an easier 
understanding of this space we suggest that reader read [3] and [4] first. 

2.MAIN RESULTS 

Before showing our main results, we need some lemmas. The following 
lemma is cited directly from [9] . 

LEMMA2. Let Q,Y, and j n be as in LEMMA1. Then, given a point 
x in Q, there is a measure a supported on Y such that 

Jp(y)da(y) = p{x) 

for all p in Vk, and 

Jd\a\ (y) < e 2«7„(fc+i). 
Now we need a famous formula. 



Stirling's Formula: n\ ^ \/2Tm(j) n . 
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Remark The approximation is very reliable even for small n. For example, 

when n = 10, the relative error is only 0.83%. The larger n is, the better the 
approximation is. For further details see [1] and [2] . 



LEMMA3. Let p x = \ and p 2 = & ^ M . Then 

V2^p k k k <kl< ^p k k k 



for all positive integer k. 

Proof. Note that ±, ^, g, ...can be expressed by±, (4) 2 , (#, (4) 4 > #) 5 > 

Now, sup 1 1, 4, 4> 4> -} = 4 im P lies that # < P2 for a11 ThusA;! ~ 

remaining part \j2~Ttp\k k < k\ follows by observing 
that V2^(l) k k k < V2^H) k ■ Vk ■ k k - fc! 



□ 



LEMMA4. Let p = ^.Then k\ < V2^p k k k ~ 1 for all k > 1. 



- vj 

e 

3 r 3 

Proo/. fc! - v / 2^(i) fe • y/k- k k = v 7 ^ • • fc*" 1 . Note that<^ *+ : fc = 1,2,3... 



3 i 3 



can be expressed by <j |, (^) 2 , (^f-) 3 > (^r) 4 "- f Our lemma follows by noting that 

C 3 13 

□ 



LEMMA5. Let fc(aj) = e - ^l x l ,/3 > 0, be the Gaussian function in i?™, 
and p be the measure defined in (4) . For any positive even integer k, 

f . . ^.h. , . n + 1 fc+n + 2 fc + n — 1 „/e . fc+n — 3 . -i . 

/ d/*(0 < -n-a n - 2~2~ ■ p—^ -02 -{k + n- l)s~ ■ (2 + ±) 



for add n, and 



o i ^it- , / ,_n n + 1 fc+n + 3 fc + n — 2 fc , . fc+n 

/ l£l < ■ n ■ a n • 2^— • • /?2 . (k + n - 2)~ 
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for even n, where p = and a n is the volume of the unit ball in R n . 



Proof. /KfdMO 

= I ^ 



fc ^n — 1 
4F 



= (|)i • n • a n • (2^) fc+ " • Jo 00 

= f vrf • n • a n • • /?§ • (*±g=2)(*±S=2 - 1) . . . (I) -^du if „ is odd, 

" 1 vrf •n-a n -2"" 1+fe -/3f(^±|^)! if n is even. 



Let b = \ k± ^ 1 ] . Then 6! < ^p h ■ b b -\ 
For odd n,b= Thus 



/ — — — fc— n+3 fc + n— 1 fc + n-3 

= v 27r -2 2 • p 2 . + n — 1) 2 . 



and 

~ . , . L> , , , , n+1 . -i . fc + n + 2 ^ fc fc+n — 1 . . fc+n — 3 

/ \Z\ k dfi(0 < vr— • n ■ a n • (2 + \) ■ 2—^- ■ f3* ■ p~^~ ■ (k + n - , 

by noting that / °° ~^du < 2 + \. 
For even n, b = k + n ~ 2 . Thus 



b\ < V^p^ • 



, — fc — n+4 fc + n — 2 , . fc+n — 4 



and 
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/, , , , x n + 1 ro-| _ K " £ . . K-\ II 

\£\ dn(£) < tt~ -n-a n - 2^2— . pi . p~^— . (k + n- 2)~ 



fc + n+3 



fc + re — 2 



fc + n — 4 



□ 



Theorem 1 Lei h(x) = e7$\ x ^ , /3 > 0, be the Gaussian function in R n , and p 
be the measure defined in (4) . Then, given a positive number bo, there are positive 
constants So, c, and C for which the following is true: If f € Ch, m and s is the h 
spline that interpolates f on a subset X of R n , then 



\f(x)-8(x)\<A"(C 6 y. 



" / N — n±± , ./o x n — A it n— i 

,where A = it 4 . (n ■ a n ) 2 • 2 4 (-^j 4 for even n and A = it 4 . (n . 

a„)5 ■ (2 + 1)3 • 2? ■ for add n, holds for all x in a cube E provided 

that (i) E has side b and b > bo, (ii) < S < 6q, and (m)every subcube of £ 
of side S contains a point of X. Here, a n deontes the volume of the unit 
ball in R n . 



The number c is equal to where 7 n is defined in LEMMAl.The 

®ln 



number C is equal to ^34 • e • y/2p/3 ■ s/n- e 2 ™ 7 ™^ • 6 3 • j n , where p 
Moreover, So can be denned by 



= vj 

e 



So = min { — 5 ,S n 

(3i-e- v ^- v / n-c 2 "T») 4 -()o 3 -7„ 



where 



& ifn = l, 



.-. t~t\ if n is odd, and n > 1, 
— - if n is even, and n > 2, 



27„(n-2) 

* if " = 2 - 



Proof. Let <5o be as in the theorem. If <5 < So, thenJ < ^-and^^ < 1. There 
integer k > 1 such that 1 < (^)fe < 2. This implies 



is an 



2 
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Now, let x be any point of the cube E and recall that Theorem4.2 of 
[8] implies that 

\f(x)-s(x)\<c k -\\f\\ h - f\y-x\ k d\a\(y) (6) 

whenever k > 0, where a is any measure supported on X such that 

J p(y)da(y) = p(x) (7) 

for all polynomials p in Vk-i- Here 



whenever k > 0. 



LEMMA5. now applies. 



For odd n, 



<* = -h \ J\t\ u Mt)\ 

KR n ) 



1 ( n+l 2k + n + 2 2fe+n-l . , „ 2fc+n-3 , i . ~l 2 

^ H 2~ s~ • p—^ .[3 k .{2k + n- 1)-^*— ■ (2+|)j 

1 f i+l , \ 1 ™+ 2 , , „ i . . 2fc+n-3 / 71 

< £ {vr— • (n • a n )2 . 2~ ■ p~ ■ (VW) fe ■ (2A+ n - l)^ - • ^2 + ± J 



n+l 1 n+2 n — 1 . , 2fc+n — 3 / ., 

* — JrJ/J " ^ ^ LEMMA 3 
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n—l , .In n—l t 2fc+n— 3 . / i . 

= vr— • (n • a„)2 • 24 • • p 3 * • (£) • ( 2 fe + n - 1)^J— • (^2 + ±) 
where /?, = ^ 

. // j 1 , , 2fc+n— 3 

= A .pf.^.Qk + n-l)—*- 

. // /T i~ n — l . . 1 n n — l 

where A = w 2 + | • 7r 4 ■ (n ■ a n )2 ■ 2* ■ p 4 . 

To obtain the desired bound on \f(x) — s(x) \ , it suffices to find a suit- 
able bound for 

1 = c k f \y~ x\ k d\a\ (y). 

R n 

Let Q be any cube which contains x, has side J n k8, and is contained 
in E. Subdivide Q into {^ n k) n congruent subcubes of side S. Since each 
of these subcubes must contain a point of X, select a point of X from 
each subcube and call the resulting discrete set Y. By Lemmal we may 
conclude that there is a measure a supported on Y which satis fies (6) 
and enjoys the estimate 

/ d\a\ (y) < e 2n ^ k . (8) 

We use this measure in (6) to obtain an estimate on /. 

Using (8) and the fact that support of a is contained in Q whose 
diameter is y/n^ n k5 we may write 

/ < A" • p k ■ £ ■ (2k + n - l)^ 1 • {^i ln k5) k ■ e 2n ^ k 

= A" • (B' ln 5) k -(2k + n- l) 2 ^ where B' = p^e" ln 

< A" • {B' ln 5) k ■ (3fc)T (if k > n - 1) 
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A" • (B' ln 53i fcl ) k 



< A 



345' 7n 5(^)4 (since A; < ^) 



A- 



a 12] 
5 7 n 3 6 ^ 3 



~i k 



where 5" = 33 B' 



A- 



_bQ_ 



(since 2^5 < A;) 



"2 I : ,,4 

where i? 6g7„ 4 54 < 1 if and only if B %^ n 5 < 1 and is guaranteed by 

5 <5 



A" 



B bhJ 



b_0 

8- 4 



A" • [C<5]^ where C = b" b 3 j n andc=^. 



Note that k > n — 1 is guaranteed by n — 1 < < 



For even n, 



< 



1 f n + 1 o 2fc + n+3 2fc+n-2 . 2fc+n-4 ~| 

i|vr— -n-a n -2 — 2 — . p — 2 — . f3 k . (2/c + n - 2) — 2 — J. 



1 f , „ 1 re + 3 ra-2 , , , , „ 2fc+n-4 ~l 

^ {vr— • (na n ) 2 ■ 2~ ■ p— ■ ■ (2fc + n - 2)~ J— j 
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n+l 1 n+3 n — 2 , , 2fc + n-4 

< 7r~j~-(na n )^-2^~-p^~-(y^^) fc -(2fc+n-2) 4 



2k+n-i 

" k (2fc+n-2) 3 , A " H=I / J H+i 21=2 

' ^ p where A = ir 4 . (na„J2 -2 * ■ p * 



A" • P 3 A 



and P3 = vm 



pi 



Now, 



I = c k j \y-x\ d\a\(y) 
R n 



2fc+n-4 



< A" • p 3 k • {2k+n -H~^~ • (Vn • 7n • k ■ 5) k • e 2 "^ 



= A" • (-B' 7n <5) fc • (2k + n - 2)^n^ where = ^^e 2 ™ 7 " 



< A" • (B'-f n 5) k ■ (3k) f (whenever k > n - 2) 



A" • (B' 7n <53f fci)* 



< A • 



34i?' 7n 5(^_)4 (fc^A 



A- 



5 7 n 5 6^^ 



ll „3 



where 5 = 3< ■ B' 



A 

< A 



_*Q_ 



Since 2^5 < k 



and B" 4 6^ 7 „ 4 <5i < 1 due to 6 < S . 



This gives 
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/ < A" B" b* ln 5 



A" • [C<5]f where C = B b 3 ~f n andc=i. 



Note that k > n — 2 is guaranteed by 5 < 5q < 



We conclude that 



|/(x) - fl (s)| < A" ■ [C*]f ||/|| fc 



whenever < 5 < 5o as stated in the theorem. 

Remark The high-level error bound for Gaussians was first put forward by 
Madych and Nelson in Theorem 3. of [9] . However their proof is incomplete. In 
their theorem the radial basis function h must satisfy the key condition 



for k > 2m, where r is a real constant and p is a positive constant. As 
pointed out by them in page 102 of [9] , r = ^ if h is a Gaussian func- 
tion. Moreover, they only treated the case = 1. However(9) holds only 
when n = 1. For n > 1, (9) should be replaced by Lemma 5 of this paper. 
Consequently their high-level error bound is essentially only suitable for 



Note that Theorem 1 can be improved if n is known in advance, es- 
pecially when n = 1 or 2. For example, suppose n = 1. Let 



/ \C\ k dK0<P h k rk 



(9) 



R 1 . 
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Then for 5 < S , 



I < A" ■ (B' ln 6) k ■ (2k)^ 
= A".(B' 7n 5) k .(2k)-2 



< A" 
= A" 



V2B' rn 5(^s) h 



A • 



A • 



B 7„ 2 <52 62 



where B" = s/2B' 



where B"^ n 2 5zb^ < 1 is guaranteed by 5 < 5q, 



A • 



B -y n b 5 



A -[CSp 



where C = B j n bo and c 



47„' 



In this case, So, C, c are the same as Madych and Nelson's results with- 
out any sacrifice. 



What's noteworthy is that in Theorem 1 the parameter S is not the 
generally used fill distance. For easy use we should transform Theorem 
1 into a statement described by the fill distance. 



Let 

d(Q, X) = sup inf \y — x\ 

y( zQ xex 

be the fill distane. Observe that every cube of side 5 contains a ball 
of radius |. Thus the subcube condition in Theoreml is satisfied when 
5 = 2d(E,X). More generally, we can easily conclude the following: 
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Corollary Suppose that h satisfies the hypotheses of the theorem, £1 is a set 
which can be expressed as the union of rotations and translations of a fixed cube of 
side bo, and X is a subset of R n . Then there are positive contants do, d ,and C 
for which the following is true: 

If / G Ch, m and s is the h spline that interpolates / on a subset X of 
R n ,then 

\f(x)-s(x)\<A" [C'dp\\f\\ h 

,where A is as in Theoreml and C = 2C, d = | with C and c defined 
in Theoreml, holds for all i in a cube ECfi provided that(i) E has side b 
and b > bo, (ii) < d < do, and (Hi) every subcube of E of side 2d contains a 
point of X. Here d denotes d(Q.,X) or d(E,X) and m = since Gaussians 
are c.p.d of order 0. 

Proof. Let do = ^ and 5 = 2d where Jo is as in the theorem. If < d < do, < 
5 < 2do = 5o, By the theorem, 

|/(x)- a (x)|<A" [C6]l.\\f\\ h 

= A" [C2d|a.||/|| h 

= A" [C'd]i.\\f\\ h 

holds for all x in 0,. 



□ 



Remark The space Ch, m probably is unfamiliar to most people. It's introduced 
by Madych and Nelson in [7] and [8] .Later Luh made characterizations for it in 
[3] and [4] . Many people think that it's defined by Gelfand and Shilov's definition 
of generalized Fourier transform, and is therefore difficult to deal with. This is 
wrong. In fact, it can be characterized by Schwartz's definition of generalized Fourier 
transform. The situation is not so bad. Moreover, many people mistake Ch, m to be 
the closure of Wu and Schaback's function space which is defined in [10] . This is 
also wrong. The two spaces have very little connection. Luh also made a clarification 
for this problem. For further details, please see [5] and [6] . 
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